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ABSTRACT 

We  propose  a  sensor  that  measures  the  number,  strength,  altitude  and  velocity  of  atmospheric  turbulence  layers. 
Recent  research  has  shown  that  pupil  plane  branch  points  contain  four  independent  and  measureable  parameters 
and  that  these  four  parameters  can  be  used  to  estimate  four  independent  turbulence  layer  parameters-number, 
strength,  altitude  and  velocity-for  each  atmospheric  turbulence  layer.  Here,  we  summarized  previous  results 
and  then  demonstrate  how  these  results  allow  for  construction  of  a  turbulence  layer  sensor. 

Keywords:  adaptive  optics,  branch  points,  atmospheric  turbulence 

1.  INTRODUCTION 

Much  of  astronomy  relies  on  images  collected  through  atmospheric  turbulence.  Turbulence  imprints  a  ran¬ 
dom  spatially  varying  phase  on  the  propagating  wave,  and  this  random  phase  is  an  error  that  degrades  image 
performance.  Adaptive  optics  is  a  long  established  means  of  compensating  for  this  error. 

Traditional  adaptive  optics  compensates  for  the  errors  by  assuming  a  phase  only  disturbance  in  the  pupil  of 
the  telescope,  i.e.  a  singe  turbulence  layer  at  the  pupil.  A  deformable  mirror  is  placed  optically  conjugate  to  the 
pupil,  and  when  the  wave  is  reflected  by  the  deformable  mirror,  the  figure  of  the  mirror  is  added  to  the  phase  of 
the  optical  wave,  thus  inverting  the  atmospheric  phase  disturbance.  The  deformable  mirror's  figure  is  calculated 
by  subdividing  the  telescope’s  pupil  into  a  regular  array  of  subapertures  and  measuring  the  wave  within  each 
subaperture  using  a  wavefront  sensor.  The  corrupted  phase  is  reconstructed  from  these  measurements  and  fed 
into  a  proportional-integral  controller  that  calculates  deformable  mirror  commands.  This  long  standing  paradigm 
has  been  recently  improved  by  both  increasing  the  number  of  deformable  mirrors,  multi- conjugate  adaptive  optics 
(MCAO)  and  an  enhanced  controller  (Predictive  Fourier  Control1)*  MCAO  improves  compensation  by  assuming  a 
“best-fit"  two  layered  structure  to  the  atmospheric  turbulence2  even  when  the  atmosphere  is  broadly  distributed. 
This  common  assumption  is  borne  out  experimentally.3  Predictive  Fourier  Control1  modifies  the  control  law  to 
account  for  a  layered  atmosphere  with  each  layer  having  a  dominant  velocity  for  a  few  tens  of  seconds.  Hence, 
measurement  of  the  three  dimensional  atmospheric  turbulence  is  of  interest  because  a  priori  knowledge  of  the 
turbulence  layers  can  be  used  to  improve  adaptive  optic  compensation. 

Of  particular  interest  to  us  is  a  peculiar  interference  effect;  after  the  wave  encounters  turbulence,  as  it 
propagates  further,  it  self-intcrferes.  Periodically,  this  self-interference  creates  zeros  in  amplitude.  These  zeros 
are  known  as  branch  points. 

Here,  we  propose  a  sensor  that  uses  branch  points  to  measure  the  atmospheric  turbulence  layers.  Branch 
points,  in  any  classic  adaptive  optics  sytstem,  in  particular  one  that  uses  a  least  mean  square  reconstructor 
only  lead  to  degradation  in  adaptive  optic  system  performance4  and  arc  regarded  strictly  as  noise.  In  our  pre¬ 
vious  work,  we  have  shown  that  branch  points  contain  a  great  deal  of  information  about  the  turbulence  that 
created  them.  In  particular,  we  demonstrated  that  the  ensemble  of  pupil-plane  branch  points  carry  four  indepen¬ 
dent  parameters.*"0  These  four  independent  parameters  can  be  used  to  estimate  four  independent  atmospheric 
turbulence  layer  parameters-number,  strength,  altitude  and  velocity.8 

In  this  work,  we  describe  a  sensor  that  estimates  those  four  parameters,  describing  how  it  can  be  built 
and  implemented  at  an  observatory.  Specifically,  we  present  the  combined  package  of  hardware  and  software 
necessary  to  build  the  sensor.  We  begin  by  first  formulating  the  problem  in  Section  2  and  and  enumerating  the 
parameters  of  importance  for  adaptive  optics  in  Section  3.  In  Section  4,  we  address  a  seemingly  major  stumbling 
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block  to  successful  implementation  of  this  sensor.  Then,  in  Section  5,  we  present  our  previous  results;  this  is  the 
foundation  upon  which  this  type  of  sesnor  can  be  built.  Section  6  presents  the  turbulence  layer  sensor,  folowing 
which  Section  7  discusses  how  this  device  was  implemented  in  our  laboratory.  Finally,  Section  8  summarizes  the 
work.  In  these  presentations,  we  demonstrate  that  branch  points,  which  heretofore  has  been  regarded  strictly  as 
npicA  in  adaptive  optics  systems,  can  be  used  to  estimate  useful  information  on  the  turbulence  layers. 

2.  THE  ZEROS  OF  AN  ANALYTIC  FIELD 

Assume  that  there  is  a  monochromatic  point  source  beyond  the  atmosphere  and  that  a  wave  horn  this  source 
propagates  through  the  atmosphere.  Atmospheric  turbulence  adds  random  optical  path  length  to  the  propagating 
wave,  and  these  path  length  differences  appear  as  random  phase  fluctuations.  As  the  wave  propagates,  the 
phase  fluctuations  cause  the  wave  to  self-interferc  as  it  propagates.  This  self  interference  leads  to  amplitude 
fluctuations.  Under  the  right  conditions,  this  self-interference  is  strong  enough  to  produce  zeros  in  amplitude. 
These  zeros  in  amplitude  are  known  to  appear  when  the  log  amplitude  variance-called  the  Rytov  parameter— of 
the  propagating  wave  exceeds  0.1.  Throughout  the  paper,  the  two  terms  ‘zeros  of  the  field’  and  ‘branch  points' 
will  be  used  interchangably. 

The  field’s  zero  in  amplitude  is  both  accompanied  by  an  equivalent  non-localized  2n  circulation  in  phase,  and 
also  each  zero  is  uniquely  associated  with  either  a  positive  or  negative  polarity.  That  is,  in  the  absence  of  noise, 
when  the  phase,  d{x,  y),  is  integrated  about  a  closed  path  about  a  single  branch  point,  the  resultant  is  always 
±2ir,4’ 10  i.e. 

f  d!0(x,y)  =  ±2ir  (1) 

J  n 

where  Jfl  •  dt  is  a  line  integral  along  the  closed  curve  fl  which  encircles  a  single  branch  point.  Note  that  the  sign 
of  the  Equation  1  is  uniquely  associated  with  cither  ±1.  For  discussional  clarity,  we  call  sign  (fa  dl  <t>{x,y))  the 
polarity  of  the  zero. 

Properties  of  the  zeros  of  the  field  are  the  foundation  for  creating  a  turbulence  layer  sensor. 

Finding  the  Zeros 

Finding  the  zeros  of  the  field  is  of  utmost  importance.  In  the  pristine  world  of  theory,  this  is  trivially  done  by 
finding  the  minima  in  the  magnitude  of  the  field.  However,  in  the  real-world,  sensors  have  finite  pixel  size,  and 
with  probability  one,  this  precludes  using  zeros  to  find  branch  points.  In  the  real-world,  in  addition  to  having 
finite  size  pixels,  they  are  also  intrinsically  noisy.  In  this  case,  it  has  been  found  that  branch  points  are  best 
found  using  Equation  l.10'11 

Creation  Pairs 

It  has  been  previously  shown0  that  zeros  appear  in  pairs  of  opposite  polarity  and  that  these  pairs  separate  as 
they  propagate.  In  addition  to  the  density  of  branch  points,  intra-crcation-pair  separation  is  key  to  the  Hpgipn 
of  a  turbulence  layer  sensor. 

3.  PARAMETERS  OF  IMPORTANCE  FOR  ADAPTIVE  OPTICS  -  WHAT  A  TULS 

MUST  MEASURE 

Atmospheric  Turbulence  and  its  Moments 

Atmospheric  turbulence  causes  phase  disturbances.  It  has  been  found  that  the  phase  disturbances  are  best 

described  by  a  structure  function  with  an  ^  power  law  in  its  power  spectrum;  this  is  known  as  a  Kolmogorov 

spectrum.12  Assuming  a  Kolmogorov  spectrum  and  using  the  Rytov  approximation,  two  common  approximation 

used  in  adaptive  optics,  it  has  been  shown12,13  that  four  moments  of  the  atmospheric  structure  junction,  for 

the  most  part,  determine  adaptive  optic  performance.  These  are  the  zeroth  altitude  moment,  the  4th  altitude 
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moment,  the  §‘b  altitude  moment,  and  the  §lh  velocity  moment.  The  functional  forms  of  these  moments  for  a 
plane  wave  are 


o\  =  0.291*o7/6/  Cl(h(z))z^ds 

(Rytov  parameter) 

(2) 

t  .L  \ 

ro  —  ^0.423fco  Jq  C'n(h(z))dsj 

(Fried'*  parameter) 

(3) 

0O  m  2.91kg  j  C*(h(z))zi/3  dz 

(the  iwtplanatic  patch) 

(4) 

fc  =  0.254*g/S  /  Cl(h(s))vi/3(z)dz 

(Creenwood  frequency) 

(5) 

\Jo 


with  ko  the  wave  number,  L  either  the  distance  to  the  guide  star  or  the  top  of  the  atmosphere,  z  the  distance 
from  the  pupil  to  that  point,  h(s)  the  altitude  at  distance  :,  and  assuming  a  single-pole  servo  filter. 

In  the  approximation  of  MCAO,  turbulence  is  layered.  Hence,  the  structure  function  is  layered.  Making 
the  same  assumption  here,  we  assume  that  there  are  NtUrb  turbulence  layers  and  these  layers  are  at  altitudes 
hi  <  hi  <•-•<  hn-  The  structure  function  is  then  given  by 

Cl(h)  =  C*(ft(r))i(h  -  hi)  +  C’(/i<*))d(A  -  *2)  +  • •  ■ •  +  CUh(z))6(h  -  hN).  (6) 


The  Parameters  of  Interest  in  Adaptive  Optica  —  ro,  ffj,  Oo  and  fc 

A  turbulence  layer  sensor  must  estimate  the  four  parameters  upon  which  ro,  a|,  8q  and  fc  depend.  These  are 
the  number  of  turbulence  layers,  and  for  each  layer,  its  velocity,  strength,  and  altitude. 

4.  CREATING  BRANCH  POINTS  ON-DEMAND 

One  potential  problem  to  practical  implementation  of  a  TuL  Sensor  is  that  for  most  observatories,  no  branch 
points  exist  in  the  optical  field.  For  an  observatory  at  an  altitude  of  3 Km  or  greater,  observing  near  zenith  at 
wavelength  of  1  pm,  and  using  the  Huffnagel-Valley  atmospheric  model,  o*  «  0.1.  This  section  demonstrates 
that  this  is  not  a  problem  because  branch  points  con  be  produced  on  demand. 

4.1  Branch  Points,  On-demand 

In  principle,  all  that  is  required  to  create  branch  points  is  to  increase  the  propagation  distance.  To  wit,  suppose 
the  actual  Rytov  parmater  is  near  zero,  but  the  required  Rytov  parameter  is  "R.  >  0.  Since  <r£  is  monotonically 
increasing  with  z,  then  there  exists  As  such  that 

Tl  =  0.29 lk70/6  /  C*(z)(z  +  Az)5'6  dz 
Jo 


In  the  layered  case,  the  A,  are  fixed  yielding 

n  m  0.291  fcj/8  (C2(h)5(fr  -  hi)(h  +  Az)s/a  +  C$(h)6(h  -  h2)(h  +  Az)s/6 . . . 

+C*(A)d(A -/»*)</» +  Az)5/a)  (7) 

the  desired  result. 
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4.2  Optically  Increasing  Propagation  Distance 


Unfortunately.  A:  may  be  kilometers,  so  direct  implementation  on  an  optical  bench  is  problematic.  However,  as 
is  well  known,  propagating  waves  of  the  same  Fresnel  number  are  diffract  ivoly  the  same.  The  Fresnel  number  of 


an  optical  wave  is  Riven  by 


(8) 


with  D  the  pupil  diameter  and  A  the  wavelength  of  the  wave. 

For  two  beams  of  different  sizes.  Doutpu,  #  two  beams  are  equivalent  as  Iour  as  and  zwput  arc 

adjusted  such  that 

j-  _  (-PoiiliHif  /»)*  _  (■0|ti;iu(/-)* 

''Jiliplll 


This  fact  is  used  by  the  optical  community  in  every  adaptive-optic  optical  design.  For  instance,  consider  a 
telescope;  as  the  wave  passes  through  a  telescope,  it  is  magnified  by  the  telescope.  Let  Doutpui  be  the  diameter 
of  the  input  aperature  and  D,„put  =  MDoutput  be  the  diameter  of  the  beam  in  the  cotide  room.  If  A;0(),  is 
required  to  obtain  T  and  is  the  maximum  distance  available  in  the  coudti  room,  then  choose  M  such 

tlrnt 

Atm  pul  =  A/  Acriu|p„|  (10) 


Since  M  is  selectable  via  the  optical  design.  A;oufpU(  of  tens  of  kilometers  can  be  equivalent  to  Ac,„,,U|  of  meters. 

Note,  since  the  altitude  of  turbulence  layers  is  arbitrary,  this  technique  varies  the  propagation  distance,  and 
hence  the  Rvtov  parameter  at  will,  and  this  allows  us  to  create  branch  points  at  will. 


5.  PREVIOUS  RESULTS  -  MEASURING  THE  PARAMETERS  OF  INTEREST 

In  previous  work,  we  have  shown  that  the  branch  points  of  the  field  contain  a  rich  amount  of  information.  In 
the  foundational  work  of  this  research  thread.9  it  was  shown  that  branch  points  arc  created  in  pairs  of  opposite 
polarity  infinitesimally  close  together.  Based  on  causality,  it  was  speculated  that  branch  points  drift  apart  as 
they  propagate,  that  the  separation  within  creation  pairs  yields  the  distance  to  the  turbulence  layer  that  created 
this  creation  pair,  and  that  the  creation  pairs  have  the  velocity  of  the  turbulence  layer  that  created  them.  This 
is  shown  graphically  in  Figttre  1.  Later  work  expanded  on  this  initial  result  greatly:  this  section  reviews  that 
work. 


"The  Creation  and  Evolution  o I  Branch  Points' 
Implications  ol  Causality 

*  Branch  points  drift  apart  as  they  propagate 

•  Branch  point  separation  is  proportional  to 
the  distance  to  the  turbulence  layer  that 
created  them 


Experimental 

Oata 

- 1  • 


CEBP  results 


o 


Creation  pairs  have  the  velocity  of  the 
turbulence  layer  that  created  them 


FiRurc  1.  .-I  gnphical  view  of  the  results  of  the  foundational  work  fortius  research  thread. 


The  Polarity  Sensor 

In  our  previous  work,  a  thought  experiment  sensor  was  assumed  that  measures  an  optical  wave  at  time  tk  and 
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returns  a  sensor  frame  (array)  of  data  that  contains  ±1  at  the  branch  point  locations  and  zeros  elsewhere.  Call 
this  the  polarity  array  nnd  denote  it  p(x,y,t).  All  the  information  needed  for  a  turbulence  layer  sensor  can  be 
extracted  from  p(x,y,t).  This  section  enumerates  the  four  independent  parameters  that  are  contained  witlun 

p(x,  (/.<)■ 


5.1  The  Velocity  Distribution 

The  velocity  distribution  is  created  from  the  set  of  all  possible  velocities  in  p(x.  yJk).  Ihe  set  of  all  possible 
velocities  is'  the  set  of  all  calculated  velocities  for  all  permutations  of  branch  point  velocities  going  forward  in 
time  for  branch  points  of  similar  polarity.5  Then  create  a  function.  D(t').  that  takes  the  elements  of  the  set  of  all 
possible  velocities  as  input  nnd  returns  the  number  of  times  each  velocity  occurs;  D(v)  is  the  velocity  distribution 
for  this  particular  velocity  set. 

FYom  D(v),  two  parameters  can  be  extracted,  the  number  of  velocity  classes  and  the  mean  velocity  of  each 
class. 


5.2  The  Velocities  Classes  in  p(x,y,t) 

Branch  points  caused  by  the  same  ntmosphcric  turbulence  layer  create  self-correlating  peaks  in  D(u).5  In  the 
ideal  case,  these  self-correlating  peaks  appear  as  delta  functions  in  D(v).  Specifically,  say  there  is  a  single 
turbulence  layer  with  Nbp  branch  points,  D(v)  will  contain  A ’bp  delta  ftmetions.  Further,  assume  that  each 
branch  point  is  measured  Br  times;  in  this  case,  the  delta  functions  will  have  n  coefficient  with  amplitude  F(BP). 
The  velocity  distribution  then  coutnins  this  term: 

•v, 

YlnBjWv-vj).  (a) 

j=t 

where  Vj  is  the  location  tliat  the  self-correlating  peak  occurs. 

Each  of  these  delta  functions  is  a  velocity  clnss  of  the  zeros  of  the  field.  In  real-world  data  (sec  Section  7), 
due  to  noise  and  cross-terms  in  the  distribution,  these  delta  functions  appear  as  large  amplitude,  sharply  spiked 
peaks. 

5.3  The  Number  of  Velocities  Classes  and  the  Velocity  of  each  Class 

The  Velocity 

The  velocities  of  the  zeroes  of  the  field  are  the  locations  of  the  delta  functions  in  D{v).  In  theoretical  formulations, 
the  velocities  are  the  v}  found  in  the  delta  functions  and  can  be  read  off  directly  from  £'vi,  T(B})S(v  —  v}).  In 
real-world  data,  the  velocities  have  to  be  estimated  from  t  he  measured  velocity  distribution.  This  is  done  using 

{Vi}  =  {locimaxu*  [Dm'a’(vr))))  /{0}  (12) 

where  the  {V<}  are  the  velocities  of  each  velocity  class,  mux\ac(-)  returns  local  maxima,  loc(max(-))  is  the  velocity 
at  which  those  maxima  occur,  and  /{0}  eliminates  the  zero  due  to  the  fixed  pattern  noise. 

The  Number 

In  theoretical  formulations,  the  number  of  velocity  classes  is  obtained  by  counting  the  number  of  delta  functions 
in  D(v).  In  real-word  data,  the  number  of  velocity  clnsscs  is  given  by 

-Vr.c.  =  dim  ({V, ))  (13) 

where  is  the  number  of  velocity  classes  and  ditn(-)  is  the  number  of  delta  functions  in  the  velocity  distribution 
(the  dimension  of  the  set  {V’,)). 

The  Importance  of  Velocity  Classes 

Finding  the  velocity  classes  is  of  fundamental  importance  to  creation  of  the  turbulence  layer  sensor.  Velocity 


classes  not  only  give  us  the  two  parameters  above,  but  just  as  importantly,  they  allow  us  to  sort  branch  points 
into  velocity  classes.  Sorting  branch  points  into  a  velocity  class  allows  us  to  measure  both  density  and  creation 
pair  separation  for  that  class. 

A  cartoon  showing  why  separation  into  velocity  classes  is  possible  is  presented  in  Figure  2.  Laboratory  data, 
in  this  case  for  velocity  filtering,  is  shown  in  Figure  3. 


(<0  (e)  (0 

Figure  2.  Cartoon  of  the  steps  in  the  identification  of  persistent  pupil  plane  branch  point  pairs.  Red  and  ■; reen  dots  indicate 
the  positive  and  negative  circulations  respectively,  (a)  The  empty  jmlanty  array,  (h)  The  polarity  arrey  with  ±1  ’s  at  the 
locations  of  positive  and  negative  circulations  respectively,  (cj  Diagram  showing  the  estimated  velocity  ■>  from  the  positive 
branch  point  locations,  (d)  Blue  lines  indicate  the  initial  pairings  following  frame  by  frame  use  of  a  walking  algorithm, 
(e)  Pairings  fallowing  the  use  of  piston  shifting  technique.  (6)  Overlay  of  measured  velocity  to  incorporate  3rd  dimension 
of  WFS  data. 


5.4  Creation  Pair  Separation  (A) 

Each  branch  point  is  placed  in  its  velocity  class.  Then,  for  each  velocity  class,  let  A,  be  tin  mean  separation 
between  creation  pairs.  It  has  been  shown8, 11  that  the  inter-pair  separation  decreases  in  a  well  defined  way  with 
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Figure  3.  .4  polarity  array  shown  in  3-D  (left)  and  a  blowup  of  the  output  of  the  pairing  algorithm  (right).  This  type  of 
processing  allows  us  to  place  h  touch  points  into  velocity  classes,  of  fundamental  importance  for  finding  -A  and  Si. 


turbulence  strength  and  propagation  distance.  For  the  the  turbulence  lnycrs  in  our  laboratory 

A,*  (Cl,VW'0z;"'l\  (14) 

with  o,  the  propagation  distance  to  turbulence  layer  i  and  the  turbulence  strength  of  layer  / 

5.5  Intra-crcation  Pair  Separation  (d) 

For  each  velocity  class,  let  <5,  be  the  mean  separation  between  the  positive  and  negative  branch  i  oints  within  a 
creation  pair.  It  has  been  shown8  that  the  intra-pair  separation  increases  in  a  well  defined  way  with  propagation 
distance.  For  the  turbulence  layers  in  our  laboratory,  we  have  shown 


d,  CX  y/Ti 


(15) 


where  the  propagation  distance  to  turbulence  layer  i.  The  proportionality  to  N/r7  matches  a  much  earlier 
observation  by  Tatarskii1'  and  gives  tis  great  confidence  in  the  veracity  of  these  results.  Our  experimental  work 
is  on-going,  hut  a  finer  estimate  at  the  current  time  appenrs  to  be 

St  *  -  :o.i)u-  (16) 

with  as  before  and  ?o.i  the  distance  to  the  onset  of  brnnch  points  for  turbulence  Inver  i. 

6.  THE  TURBULENCE  LAYER  SENSOR 

Conceptually,  this  seusor  is  simple.  It  takes  ns  input  an  optictil  wave  in  t he  pupil  plane  of  a  telescope  and 
returns  as  output  the  number,  velocity,  altitude  and  strength  of  the  atmospheric  turbulence  lay.-rs.  Since  four 
independent  turbulence  layer  parameters  arc  desired  to  be  found,  four  independent  parameters,  related  to  the 
output  parameters,  must  be  found.  (Note  that  the  four  parameters  contain,  in  reality.  4ArP.r.  ~  1  variables.) 
Specifically. 


Turbulence  Parameters 
#  of  turbulence  layers  Aiurt 

velocity  of  the  turbulence  layers  V, 
altitude  of  the  turbulence  layers  h, 
strength  of  the  turbulence  layers  C„., 


Pupil-plane  Parameters 
•V.  c  #  of  velocity  classes 

v,  mean  of  each  velocity  clas.- 

A,  branch  point  density 

St  intra-pair  separation 


Whence  the  four  pupil  plance  parameters  was  described  in  Section  5.  The  purpose  of  this  section  is  to  describe 
the  device  that  measures  the  four  pupil  plance  iniramelers,  mid  how  that  device  converts  those  parameters 
into  useful  information.  This  includes  a  description  of  how  the  sensor  would  be  operationally  used,  the  optics 
necessary,  the  algorithms  that  convert  the  raw  wavefront  sensor  data  into  the  polarity  array,  the  algorithms  that 
translate  the  polarity  into  meaningful  output  data,  and  TuLS  at  your  observatory. 


6.1  Cartoon  and  Description  of  The  Sensor’s  Components 

The  turbulence  layer  sensor  can  be  broken  into  tltreo  major  components,  the  optics,  t he  polarit .  sensor,  and 
the  processing  algorithms.  These  arc  shown  pictoriallv  in  Figure  4.  The  optics  capture  the  input  light,  add 
propagation  distance  as  required,  and  relay  it  to  the  polarity  sensor.  The  polarity  sensor  is  comprised  of  a 
wavefront  sensor  with  additional  processing;  its  output  is  the  polarity  array  at  each  propagation  distance.  The 
polarity  arrays  arc  passed  to  the  processing  system  where  the  turbulence  layer  parameters  are  cal.  ulated. 


Figure  4.  .4  four  block  diagram  of  the  TuL  sensor  (shown  in  gray). 


The  sensor  takes  time  resolved  measurements  of  the  field  at  different  propagation  distances.  The  purpose  of 
the  measurements  is  to  take  enough  duta  to  sort  branch  poitits  into  velocity  classes  and  measure  A,  and  6,  for 
cncli  velocity  class.  To  do  this,  firstly,  the  sensor  creates  branch  points.  At  the  cartoon  level.  Figure  5  shows 
pictoriallv  how  this  is  done.  A  point  source  is  shown  above  the  atmosphere,  propagating  through  it  to  the  pupil 
of  the  telescope.  The  TuLS  is  shown  in  gray.  It  takes  the  input  wave  and  increases  the  propagation  distance  until 
branch  points  form;  the  branch  points  arc  shown  ns  blue  and  green  circles.  The  propagation  distance  is  then 
repeatedly  increased  (shown  jis  the  ;,)  with  time-resolved  measurements  of  the  field  taken  at  each  propagation 
distance. 


0  •  I. 

0  o  I, 

TiA  Motor 

Figure  5.  A  pictorial  view  of  how  the  turbulence  layer  sensor  varies  its  parameter  spare;  it  repeated  increases  propagation 
distance,  taking  time-resoved  measurements  of  the  field  at  each  distance.  The  '1\iLS  is  shown  as  the  gray  tux  beneath  the 
telescope 's  pupil. 


6.2  The  Optics 

Other  than  simple  relays,  the  optics  must  accomplish  two  tasks,  first,  it  must  cookie-cut  the  pup.l  to  a  proper 
dimension  (see  below)  and  second,  it  must  create  addition!  propagation  ns  required  by  the  alorithm-.  To  do  this, 
the  optics  first  capture  the  input  beam  passed  by  the  primary  optical  beam  train.  For  simplicity,  will  assume 
that  a  pupil  is  relayed  to  the  TuLS  input  aperture.  (This  is  not  really  required  for  the  TuLS  to  work;  however,  if 
this  is  not  the  case,  additional  optics  are  required,  trivial  to  implement  but  not  addressed  here.)  A  piece  of  the 
input  aperture  is  selected  to  be  interrogated.  This  piece  of  the  pupil  is  relayed  to  an  optical  trombone  which  adds 
additional  propagation  jus  required  by  the  algorithms.  It  is  then  relayed  to  a  wavefront  sensor  which  measures 
the  field  and  outputs  the  sampled  wavefront  at  known  time  intervals  sufficiently  small  in  order  time  resolve 
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the  field.  This  is  shown  pictorially  in  Figure  (>.  The  renininder  of  the  system  is  processing  to  b.  discussed  in 
Section  6.4. 


Figure  6.  .4  cartoon  of  TuLS  showing  additional  detail  on  the  type  of  optics  nt/itired. 


Discussion  Concerning  The  Pupil 

The  TuLS  takes  as  input  a  pupil,  that  is,  a  telescope  pupil  is  relayed  to  the  input  aperture  of  the  TuLS.  For  large 
telescopes,  only  a  piece  of  the  pupil  is  necessary.  For  instance,  in  our  laboratory,  we  mock  up  n  l.'.m  pupil.  So, 
for  large  telescopes,  the  first  thing  TuLS  does  is  limit  the  aperture  to  approximately  1.0m. 

The  size  of  the  aperture  necessary  is  determined  by  the  wind  speed.  The  slower  the  wind  speed,  the  smaller 
the  interrogated  piece  of  the  aperture  must  be.  This  is  because  the  wind  speed  and  aperture  size  'let ermine  the 
number  of  times  each  branch  point  pair  is  measurcsl  prior  to  exiting  the  pupil."’ 

Discussion  of  The  Trombone 

The  trombone  is  the  optical  component  that  adds  additional  propagation  distance.  The  distauo  to  be  added 
depends  on  the  altitude  and  strength  of  the  turbulence  layers.  Since  tin*  structure  function  is  i  function  of 
altitude,  given  by  an  atmospheric  model,  for  instance  the  Huffnagel*  Valley  model,  the  strength  dependence  can 
be  obtained  once  the  altitude  is  kuown.  In  output  space,  for  our  laboratory,  all  additional  propagation  was 
less  than  20 A" m  and  for  most  runs,  it  was  less  then  12A*m.  Given  these  distances,  once  an  input  beam  size  is 
optically  fixed,  the  length  of  the  trombone  is  determined  by  Equation  9.  For  our  benchtop  implementation,  the 
trombone  moves  less  than  140cm. 

The  presicion  to  which  the  ou-bench  trombone  must  be  moved  is  a  function  of  Equations  9  and  10.  For  our 
benchtop  implementation,  trombone  movement  of  1mm  is  sufficient. 

Discussion  of  The  Polarity  Sensor 

The  polarity  sensor  is  composed  a  wavefront  sensor  with  additional  processing  that  calculate'  the  polarity 
array.  The  wavefront  sensor  can  be  any  standard  design,  for  instance,  a  Shack-Hart  maim  or  >•  If- referencing 
interferometer.  It  can  operate  iu  any  band  from  the  visible  to  the  near  IR.  but  operation  in  ti  e  visible  will 
decrease  the  trombone  length.  The  number  of  wavefront  sensor  subnpert tires  depends  on  the  pupil  size;  for  our 
laboratory  data,  a  sensor  with  200x200  pixels  was  used  on  an  output  aperture  size  of  1.5m.  The  frame  rntc  of  the 
camera  is  determined  by  the  wind  velocity;  take  enough  samples  per  the  wind's  rorellatiou  time  to  resolve  the 
motion.  Note,  wavefront  sensors  currently  installed  at  astronomical  observatories  will  have  difficulty  detecting 
the  branch  points  since  the  supnperturc  size  is  so  large. 

6.3  The  Parameter  Sets  Required  to  Estimate  the  Turbulence  Parameters 

The  TuLS  requires  time  resolved  open  loop  samples  of  the  field.  Tin*  number  of  samples  depends  on  the  depends 
on  the  wind  speed,  the  size  of  the  pupil  mapped  to  the  WFS,  ns  well  as,  noise.  For  our  laboratory  data,  200 


9 


uncorrelated  measurements  at  each  propngut ion  distance  was  sufficient.  Data  measured  at  tea  propagation 
distances  after  the  onset  of  branch  points  is  sufiicient  to  get  good  estimates. 

6.4  Processing 

There  are  six  indpeudent  algorithms  that  must  be  run  to  calculate  the  four  pupil  plane  parameter.  These  will 
be  discussed  in  turn. 

Assume  that  the  wavefront  sensor  has  measured  a  sufficient  number  of  frames  of  dnta  at  cad.  propagation 
distance.  Then,  for  each  propagation  distance,  the  following  algorithms  are  run  in  the  order  given.  References  to 
the  previous  paper  that  describes  the  results  are  annotated.  The  first  algorithm  finds  the  location  and  polarity 
of  the  branch  points,  loads  the  polarity  array,  and  calculates  the  overall  branch  point  density  (si  e5,6-10).  The 
second  algorithm  calculates  the  velocity  distribution  (see').  The  third  algorithm  establishes  an  initial  branch 
point  pairing  using  in-frame  branch  cuts  (see7).  The  forth  algorithm  identifies  the  creation  pairs  using  the  results 
of  the  third  algorithm  and  velocity  filtering  (see’*).  The  fifth  algorithm  finds  the  velocity  classes  of  the  polarity 
array  (see16).  The  sixth  algorithm  equates  a  velocity  with  each  unique  subset  of  p(x.y,t)  ami  the  creation 
pairs  are  assigned  to  a  unique  subset  (see7'’').  Once  these  six  algorithms  have  completed,  we  haw  in-hand  the 
polarity  array,  and  the  velocity  distribution  for  this  polarity  array.  As  described  in  Section  5.  from  these,  four 
independent  pupil  plane  parameters  can  he  calculated.  These  are 

1.  the  number  of  velocity  classes  contained  in  the  velocity  distribution, 

2.  the  mean  velocity  of  each  velocity  class. 

3.  the  creation  pair  separation  for  each  velocity  class,  and 

4.  the  intra-pair  separation  for  each  velocity  class. 

As  is  stated  in.6  these  four  parameters  are  all  that  is  required  to  calculate  four  turbulence  layer  parameters. 
Turbulence  Layer  Estimation  Number 

The  number  of  turbulence  layers,  N,urk  is  equal  to  the  number  of  velocity  classes,  .V,,  c.. 

Afurfc  =  A  i  (17) 


Turbulence  Layer  Estimation  —  Velocity 

The  velocity  of  each  turbulence  layer.  V’,,  is  the  scaled  velocity  of  each  velocity  class  in  p[x,y,  t). 

V,  =  C„  <-,  (18) 

where  i  denotes  the  i*h  layer,  n,  is  the  measured  speed  in  pixels/frame,  V,  is  the  turbulence  velocity  in  me¬ 
ters/second  and  Cv  is  a  constant  particular  to  each  experimental  setup.  For  our  testbed,  Cv  =  100ft  ( 1A)  =  6.25. 


Turbulence  Layer  Estimation  Altitude  unit  St  length 

The  strength  and  altitude  of  each  turbulence  layer  is  found  by  lit  verting  the  A  and  S  curves  for  each  velocity 
class. 

For  data  taken  in  our  laboratory.*  the  turbulence  strength  is  given  ns 


r’i  _ [it/J  j a 

t-n.i  1  ”1  “i 


where  c  is  a  constant,  and  the  altitude  of  the  layer  is 


(19) 


;,=  !+</  <$r‘  (C* ,) 


!»/’» 


(20) 


where  d  is  a  constant  and  Co  the  distance  to  the  onset  of  branch  point  formation.  We  make  no  claim,  as  yet, 
that  these  results  are  definitive,  merely  that  they  are  suggestive  that  well-defined  curves  exists  in  nature,  and 
these  curves  can  be  inverted. 
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6.5  TuLS  at  your  Observatory 

While  the  previous  .sections  outlined  theoretical  considerations  for  the  TiiLS,  this  section  covers  practical  issues. 
HTitrt  in  the  system  would  it  (could  it)  go) 

The  TuLS  requires  open  loop  data.  Hence,  a  logical  place  would  be  prior  to  the  steering  mirro:  (we  assume 
that  standard  adaptive  optics  system:  the  beam  encounters  a  steering  mirror  followed  bv  the  dcfoi  mable  mirror 
followed  hv  the  wavefront  sensor).  Also,  this  device  would  best  be  permanently  installed  on  the  optical  system: 
this  is  not  a  requirement,  hut  would  reduce  down  time  due  to  realignment  at  each  data  collection. 

How  does  it  get  light  (pick  off  considerations )'? 

A  slide-bnsed  pickofF mirror  could  be  placed  in  the  optical  system  and  inserted  for  TuLS  data  collection;  a  dichroic 
pickoflf  would  also  he  useable  here.  However,  it  could  also  bo  placed  after  the  deformable  mirror,  as  long  ns,  the 
system  is  running  open  loop  while  TuLS  is  collecting  data. 

Isolation  requirements 

A  TuLS  sensor  would  have  the  same  requirements  ns  a  wavefront  sensor,  i.e.  a  vibration  stabiliz'  d  part  of  the 
l>enin  train.  Temperature  fluctuations  would  not  be  a  consideration  outside  of  this,  that  is,  change-  to  trombone 
length  would  not  affect  successful  operation. 

How  big  is  it) 

A  design  could  ettsilv  fit  TuLS  within  a  2mx2 hi  part  of  the  optical  bench.  This  could  be  reduced  it  bench  space 
is  at  a  preimium. 

H7itn  to  use  it) 

The  device  would  be  used  often  enough  to  determine  bulk  wind  velocity.  This  would  he  a  few  times  per  night 
in  the  gcuend  vicinity  of  the  science  object.  Also,  the  data  collection  would  take  no  more  than  a  few  minutes. 
At  present,  the  processing  takes  a  few  hours  using  off-the-shelf  PCs:  investment  in  dedicated  pn-essors  could 
shorteu  this  to  minutes. 

The  Sensor  in  Operation  How  do  you  use  it) 

To  be  used,  first  a  bright  unresolved  star  must  be  found  in  the  general  vicinity  of  the  science  ohjc  t.  Insert  the 
flip  mirror,  redirecting  the  light  to  TuLS.  Do  a  quick  scan  for  turbulence  producing  layers,  finding  the  general 
altitudes  of  the  strong  layers.  Set  TuLS  at  the  altitude  (A;, of  the  lowest  turbulence  layer,  then  take  data 
in  increasing  .r>00m  increments  from  that  point.  Finally,  run  the  analysis  scripts.  The  total  time  i  •  do  this  will 
be  determined  by  the  bright  ness  of  the  source  aud  the  speed  of  the  turbulence,  hut  it  should  he  n.  more  than  a 
few  minutes. 

7.  IMPLEMENTATION  OF  TULS  IN  ASALT  AND  RESULTS  FOR  THAT 

IMPLEMENTATION 

A  version  of  the  TuL  sensor  was  implemented  in  the  ASALT  laboratory.  This  section  presents  how  1  his  was  done 
along  with  results  from  that  implementation. 

7.1  The  Optics  and  Optical  Hardware 

Instrumentation 

The  data  presented  here  is  taken  on  a  testbed  equipped  with  a  well  calibrated  atmospheric  turhulcn.  >•  generator17 
consisting  of  two  G-inch  disks,  each  etched  with  a  Kolmogorov  structure  function.  To  this  standard  -imulator,  we 
tiddcd  a  variable  length  opticid  trombone.  These  are  shown  in  Figure  7.  There  we  see.  in  the  low.  r  left  corner, 
the  turbulence  simulator,  followed  by  the  trombone’s  pickoff  mirror,  followed  immediately  by  the  1 1  unbone.  The 
next,  the  standard  components  of  an  adaptive  optics  sytem  are  seen,  the  steering  mirror,  the  deformable  mirror 
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and  the  wavefront  sensor.  Due  to  our  limited  budget,  the  adaptive  optic's  wavefront  sensor  was  u-.-d  in  lieu  of  a 
dedicated  wavefront  sensor:  this  is  possible  because  our  standard  wavefront  sensor  is  eight  times  « .versampled. 

A  typical  parameter  set  for  a  velocity  measurement  is  shown  in  Table  1.  For  purposes  here,  only  t  wo  items  are 
of  note.  First,  in  column  two  the  number  of  frames  is  200,  that  is.  only  two  hundred  independent  measurements 
are  required  for  us  to  obtain  estimate  the  turbulence,  in  real-world  data,  this  would  be  tens  of  seconds  since  the 
branch  points  must  cross  wavefront  sensor  suhapertures  in  order  to  he  observed.  Second,  those  tests  have  two 
turbulence  layers.  Also  (obviously),  for  full  turbulence  layer  estimation,  variation  in  propagation  distance  would 
also  have  to  be  done. 
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Table  1.  Tiirbulencc  parameters  used  to  demonstrate  both  that  our  method  is  invariant  with  reaped  to  Turl’nlence  strength 
and  also  that  the  assumption  of  branch  point  persistence  has  experimental  basis.  Bold  print  is  used  to  identify  the  branch 
point  producing  layers. 


7.2  Experimental  Results 

For  the  parameters  in  Table  1.  datn  was  collected  and  the  polnritv  army  calculated.  Projections  ..f  the  polarity 
arrays  onto  the  a ■  -  I  plane  arc  plotted  in  Figs.  8  plots  (a),  (c),  and  (e):  branch  point  trajectories  through  the 
array  are  readily  apparent.  The  velocity  distributions,  D{vx)  arc  shown  in  Figs.  8  plots  (b).  (d).  and  (f).  The 
peaks,  characteristic  of  turbulence  layers,  are  readily  apparent. 

Number  and  Velocity  of  Turbulence  Layers 

Focusing  on  the  delta  functions  in  the  velocity  distributions,  we  see  that  the  number  of  turhuleno  layers  is  one, 
one  and  two.  for  Configurations  1.  2.  and  3.  respectively,  just  its  was  set  in  Table  1.  Turning  our  attention  to 
the  location  of  the  delta  function,  for  Configuration  1.  the  single  strong  peak  is  centered  at  0.5  pixel/ frame, 
for  Configuration  2.  the  single  strong  peak  is  centered  at  —1.5  pixel/ frame,  for  Configuration  3.  tiio  two  strong 
peaks  are  centered  at  0.5.  —1.5  pixel/ frame,  just  ns  was  set  in  Table  1. 

Strength  and  Distance  of  the  Turbulence  Layers 

Many  other  datasets  have  been  taken;  Table  2  tahularizcs  the  parameter  sets  for  one  such  data  run.  In  these 
runs  the  turbulence  was  set  at  eleven  strengths,  and  at  each  of  these  strengths  the  propagation  of  the  layer 
to  the  measurement  device  was  increased.  Results  for  this  representative  datn  set  showing  t In*  variation  in 
propagation-necessary  to  estimate  A  and  <5— are  shown  in  Figure  0.  Each  filled  circle  represent.-  ;»  data  point. 
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Figure  7.  The  turbulence  layer  sensor  implemented  on  our  testbed,  the  opto- mechanical  CAD  drau-iim  of  the  optical 
hardware  as  it  appears  on  our  bench. 
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For  these  datasets  and  using  the  techniques  already  mentioned,  the  curves  were  inverted  to  yield  Equation  19 
for  turbulence  strength  and  Equation20  for  turbulence  altitude. 
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data  was  collected  for  each  of  11  positions  of  the  optical  trombone,  ranging  in  5  cm  increments  from  0  t-  50  cm.  This 
corresponds  to  ranges  of  propagation  distances  and  Hytov  parameters  for  each  Ensemble.  The  term  ’pun. '  refers  to  the 
clear  center  portion  of  the  phase  screen.  To  be  placed  at  ’puck’  means  the  light  was  passing  through  the  cent t  r  of  the  screen 
and  not  contributing  to  the  turbulence  simulated  by  the  ATS. 
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Figure  8.  Initial  demonstration  data  for  Table  l.  lejt  column:  x-t  projections  of  p(x.y.t).  Right  column:  matching 
velocity  distributions,  D(ur).  Top  row:  Table  Configuration  I:  r0  =  6.91,  =0.19.  Middle  row :  Tabic  C'-nfiguration  3: 

rO  =  8.12,  o-j  =  0.29.  Bottom  row:  Table  Configumtion  3:  rO  =  6.25,  <7j  =  0.-1”.  In  each  case  both  phase  wheels 
were  spinning,  iVnc h  =  (1,-3).  In  Config.  I  the  branch  points  are  associated  with  the  low  altitude  phase  wheel  and 
the  spike  is  at  0.5 pixels/ frame.  In  Config.  3  they  arc  associated  with  the  high  altitude  phase  wheel  and  the  spikes  is  at 
—l.apixels/  frame.  In  Configuration  3,  branch  points  form  from  both  layers  and  two  non-zeros  spikes  at-  seen  at  —1.5 
and  0.5 pixels/ frame  respectively. 


(»)  (b) 

Figure  9.  Density  is  shown  in  plot  (a),  and  separation  in  plot  (b).  Both  are  scaled  and  normalized  j-'operly. 
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8.  DISCUSSION  AND  SUMMARY 

Here,  \vc  proposed  a  sensor  that  measures  the  number,  strength,  altitude  and  velocity  of  atmospheric  turbulence 
layers.  This  is  based  on  our  previous  work  linked  the  zeros  of  the  field  to  the  turbulence  layer  that  •  rented  them. 
Our  other  recent  research  has  shown  that  the  zeros  of  the  optical  wave  contain  four  independent  pa  atneters  that 
can  be  measured.  These  four  independent  parameters  can  be  used  to  estimate  four  independent  tmhulcncc  layer 
parameters— number,  strength,  altitude  and  vclocity-for  ench  atmospheric  turbulence  layer.  Though  branch 
points  arc  not  commonly  present  in  the  observing  geometry  typically  used  in  astronomy,  it  is  trivial  to  show  that 
by  artificially  increasing  tlx*  propagation  distance  using  an  optical  trombone,  branch  points  can  i>c  created  on 
demand. 

Here  we  presented  the  rationale  for  a  sensor  that  creates  branch  points,  measures  four  independent  pupil  plane 
parameters  and  estimates  four  atmospheric  turbulence  Inyer  parameters.  We  also  presented  an  implementation 
of  this  sensor  in  the  ASALT  laboratory.  Within  the  assumptions  of  a  layered  atmosphere  and  Tavlor  flow,  the 
proposal  to  for  this  scusor  is  well  founded;  however,  more  work  is  necessary  and  on-going.  First,  these  results 
must  be  extended  to  study  the  effect  of  broadband  light  on  the  measurements.  Secondly,  the  veracity  of  these 
results  in  the  presence  of  a  distributed  atmosphere  must  be  confirmed.  Finally,  these  results  have  t..  be  validated 
in  field  experiments;  a  plan  to  do  so  hits  been  conceived,  but  a  field  device  has  not  yet  been  built.  1'his  research 
is  on-going  anti  results  will  be  presented  as  they  become  nvnilable;  in  particular,  we  invite  you  to  <>ur  upcoming 
presentation  at  IEEE  in  Dig  Sky,  Montana  where  we  will  present  wave  optic  simulation  of  this  same  effect  and 
our  related  work  in  cnutroling  a  nwd(\)  device 
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